FRAMES AND OVERSAMPLING FORMULAS 
FOR BAND LIMITED FUNCTIONS 



VINCENZA DEL PRETE 

Abstract. In this article we obtain families of frames for the space of 
functions with band in [—0, 0] by using the theory of shift-invariant spaces. 
Our results are based on the Gramian analysis of A. Ron and Z. Shen and a 
variant, due to Bownik, of their characterization of families of functions whose 
shifts form frames or Riesz bases. We give necessary and sufficient conditions 
for the translates of a finite number of functions (generators) to be a frame 
or a Riesz basis for . We also give explicit formulas for the dual generators 
and we apply them to Hilbert transform sampling and derivative sampling. 
Finally, we provide numerical experiments which support the theory. 



1. Introduction 

In many signal and image processing applications, images and signals are as- 
sumed to be band limited. A band limited signal is a function which belongs to the 
space of functions in L^(E) whose Fourier transforms have support in [—0,0]. 
Functions belonging to this space can be represented by the Whittaker-Kotelnikov- 
Shannon series, which is the expansion in terms of the orthonormal basis of trans- 
lates of the sine function. The coefficients of the expansion are the samples of the 
function at a uniform grid on M, with "density" ^. This sampling density is usually 
called the Nyquist density. 

The theory has been extended in many directions by several authors. In one of 
these extensions the sine orthonormal basis has been replaced by Riesz bases (see 
for example |Hij by J. R. Higgins) or frames. Frames generally are overcomplete 
and their expansion coefficients are not unique. Their redundancy is useful in 
applications because the reconstruction is more stable with respect to errors in the 
calculation of the coefficients and it allows the recovery of missing samples [F]. 
In signal analysis frames can be viewed as an "oversampling" with respect to the 
Nyquist density. 

The second extension consist in using more than one function to generate the space 
. In this case the Riesz basis or the frame are formed by the translates of a finite 
family of functions and the expansion formula is called a multi-channel sampling 
formula |P] pT] , 

Finally we mention that there is a huge literature on the problem of reconstruc- 
tion of signals from non-uniform samples. Since this paper is only concerned with 
uniform sampling, we refer the reader to the recent article of A. Aldroubi and 
K. Grochenig fAG| and the references given there for an extensive review of the 
problem of non-uniform sampling. 



Key words and phrases, frame, Riesz basis, shift-invariant space, sampling formulas, band 
limited functions. 
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In this paper we construct multi-channel uniform sampling formulas for band- 
limited functions using the theory of frames for shift-invariant spaces. A to-shiit- 
invariant space is a subspace of L^(IR) that is invariant under all translations 
TktaT^ € by integer multiples of a positive number to- We recall that is 
to-shift-invariant for any to- A set $ in a io-shift-invariant space S is called a set 
of generators if S is the closure of the space generated by the family E^^t^ = 
{Tkta fjip^^jk^ Z}. The space S is said to be finitely generated if it has a finite 
set of generators. Finitely generated shift-invariant spaces can have different sets 
of generators; the smallest number of generators is called the length of the space. 
The structure of finitely generated shift-invariant spaces was investigated by C. de 
Boor, R. DeVore and A. Ron with the use of fiberization techniques based on the 
range function [BDR . These authors gave conditions under which a finitely gener- 
ated shift-invariant space has a generating set satisfying properties like stability and 
orthogonality. Successively, A. Ron and Z. Shen introduced the Gramian analysis 
and extended the results of [BDR] to countable generated SI spaces [RS . In their 
paper they characterized sets of generators whose translates form Bessel sequences, 
frames and Riesz bases. For finitely generated spaces these conditions are expressed 
in terms of the eigenvalues of the Gramian matrix. In concrete cases it would be 
useful to have more explicit conditions expressed in terms of the generators or their 
Fourier transforms. In this paper, using a result of M. Bownik, we obtain these 
more explicit conditions for the space B^^ j^. We also give explicit formulas for 
the Fourier transforms of the dual generators. 

In the last part of the paper we use these results to obtain multichannel over- 
sampling formulas for band limited signals and we apply them to the derivative 
sampling. Finally, we support the theory with some numerical experiments. 
The paper is organized as follows. In Section 2 we collect some results on frames 
for shift-invariant spaces and find the range of B^ as a to-shift-invariant space. 
In Section 3 we consider a family <f> = {ipi^ ip2-, ■ ■ ■ ^n} of TV generators, where N 
is the length of B^ as a to-shift-invariant space. We give a necessary and sufficient 
condition for to be a frame or a Riesz basis for B^^. The condition is expressed 

in terms of the pre-Gramian. To prove this condition we use a result of Bownik 
which characterizes the system of translates i?$,to as being a frame or a Riesz family 
in terms of "fibers" [B]. 

In Section 4 we find the family (f>* of dual generators. Here we use the fact that, in 
the fibered representation of B^^, the frame operator is unitarily equivalent to the 
operator of multiplication by the dual Gramian matrix G'$^t^(a;) — ^ {x)J^ ^ {x) 
acting on the fiber over x. The problem of finding the dual generators is reduced 
to solving the matricial equation ^ {x) — ^ {x)J^ ^ (x) J^* ^ (x) in the un- 
known J^t J (x) in [0,h]. Where ^ (x) is invertible, the solution is the inverse 
of j {x), elsewhere it is given by the Moore-Penrose inverse of ^ (x). We give 
the expressions of the Fourier transforms of the dual generators as cross products 
of translates of the vector <I> = (p2, ■ ■ ■ (Pn)- 

In Section 5 we apply the results of the previous sections to obtain two and three- 
channel sampling formulas for functions of B,^ and we apply them to derivative 
sampling. We also provide some numerical experiments. 
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2. Preliminaries 

In this section we collect some results on frames for shift-invariant spaces to be 
used later. We begin by introducing some notation. The Fourier transform of a 
function / in L^{M.) is 



The convolution of two functions / and g is 

/ * gi^) = / fi^ - y)g{v) dy, 



/27r 



so that !F{f * g) = TJJ-g. Let /i be a positive real number; is the space of 
/i-periodic functions on K such that 



< CXD. 



(2.1) \\f\\Ll = {\f\f{xWdxy 

With the symbol £■^(2; C^) we shall denote the space of C^-valued sequences c : 
(c(n))z such that 

1/2 



l|c||£. = (El^WI' 



< oo. 



Let _ff be a subspace of L^(IR). Given a subset (f> = {fjij = 1, ■ • ■ , of H and a 
positive number to denote by E,^,to the set 

^*,to = WntoVj: nel, j = 1, . . • ,iV}. 

Here Taf{x) — f{x + a). The closure of the space generated by i?$,to will be denoted 
by S<p^ta- The family E<s,,t„ is a frame for if the operator ^ : li?(TL\ C^) H 
defined by 

N 

j=l n£Z 

is continuous, surjective and ran(T^ ^ ) is closed. The family E^^ta is a frame for 
H if and only if there exist two constants Q < A < B such that 

N 

(2.2) ^ll/lp < 5] 5] l(/,r„,„^,)P < BWff V/ e ff. 

j = l nSZ 

The constants A and B are called frame bounds, li A = B the frame is called tight 
and if A = B = 1 a Parseval frame. Denote by T| ^ : H ^ ^^(^; C^) the adjoint 
of defined by 

(2.3) (T|,,J),(n) = (/,r„t„(^,) n ^ j = I . . . , N . 

The operator ^ f. : H H is called frame operator. The set i?$,t„ is a frame 
for H if and only if the frame operator is continuously invertible and 

N 

j=l nez 



Observe that (2.2 1 can be written AI < ^ ^ < BI, where / is the identity 
operator on H. Denote by <&* the family $* = {ifij^j = 1, . . . , N}, where 

(2.4) V* = (r<,_,^T* , J- V, 1 < J < iv. 
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If E^^ta is a frame for H then , ta is also a frame (the dual frame) , and j T^, ^ = 
'^$*,to^i,t„ = ^- Exphcitly 

JV TV 
(2-5) f ^^^{f-,Tnt,^])Tnt^^] ^^^{f,Tnt,^j)Tnt^^] ^feH. 

The elements of $ are called generators and the elements of <&* rfua/ generators. If 
the family E^^ta is a frame for and the operator ^ is injective, then Et^^t^ is 
called a Riesz basis. 

In what follows to is a positive parameter. To simplify notation, throughout the 
paper we shall set 

to 

A subspace S of L^(M) is io-shift-invariant if it is invariant under all translations 
by a multiple of to. The following bracket product plays an important role in Ron 
and Shen's analysis of shift-invariant spaces. For / and g £ L'^{M.), define 

(2.6) [f,9] = hY,f{-+jh)9{-+.]h). 

Note that [f,g] is in L\ and ||[/, = II /Hi- The Fourier coefficients of [f,g] are 
given by 

(2.7) [f,9ni)^{f,m^g) iez. 

Indeed 

[Lgfie)^ £Y.^Mfm^)e-^"'^dx^ I f{xyg{x) e-^'^^^dx 

If S is a to-shift-invariant space and there exists a finite family <i> such that S — 
S^,t^ , then we say that 5 is finitely generated. Riesz bases for finitely generated 
shift- invariant spaces have been studied by various authors. In [BDRJ the authors 
give a characterization of such bases. A characterization of frames and tight frames 
also for countable sets $ has been given by Ron and Shen in [RSI]. The principal 
notions of their theory are the pre-Gramian, the Gramian and the dual Gramian 
matrices. 

The pre- Gramian ^ is the /i-periodic function mapping M to the space of oo x A^- 
matrices defined on [0, h] by 

(2.8) {J^^^^)^^{x) = ^/h^,{x + jh), j e Z, £ = I, . . . , N. 

The pre-Gramian ^ should not be confused with the matrix-valued function 
whose entries are y/h'(pi{x+jh), for all a; G M, which is not periodic. The spectrum 
of the space S^,t^ is defined as 

(2.9) a{S^,tJ = {xeR:,h^t^ix)^0} 

or, equivalently, as the support of X^jLiI^i Of course, since the functions (pj 
are defined only up to a null-set, the support is intended in the sense of distributions, 
i.e. as the complement of the largest open set on which the function X^jLiI^j^] 
vanishes as distribution. It was proved in |BDRj that the spectrum of a finitely 
generated space depends only on the space itself and not on the particular selection 
of its generators. 
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Denote by J|, ^ the adjoint of ^ . The Gramian matrix Gq,^ta — t t 
the N X N matrix whose elements are the ^.-periodic functions 

(2.10) (G$,t„),,= [^,^^]- 

The dual Gramian G^^to = J<s> t t the infinite matrix whose elements are 

AT 

(2.11) {G^^to)jg^ h'^Tjh^p^i Tihip^, jJeZ. 

n=l 

The importance of these two matrices lies in the fact that the Gramian matrix 
represents the operator T<J ^ ^ and the dual Gramian represents the operator 
T<j ( ( and many properties of these operators can be studied by looking at 
them. Indeed, by the theory of Ron and Shen [RS , after conjugating with an 
isometry, the operator ^ ^ can be decomposed into a measurable field of 
operators, acting on £^(Z), which are represented b y th e dual Gramian matrix 



G$,t^ in the canonical basis (see formulas (3.4| and (3.5l below). Similarly, after 
conjugation with a Fourier transform, the operator T| ^ ^ is represented by the 
Gramian G^^t^. 

Denote by L^(M;^^(Z)) the Hilbert space of i!^(Z)-valued /i-periodic functions on 
K such that 

(2.12) WfhimHz)) ={ll 11/(^)11^^ dx) " < oo. 

For every / in L^(IR) we denote by Chf the ^^(Z)-valued function defined on [0, h] 
by 

(2.13) Chf{x) = VhJ2fix + ^h)6t x(z[0,h] 

and extended to M as a periodic function of period h. Here {Si : £ E Z} is the canon- 
ical basis of The map / CiJ is an isometry of L'^{R) onto Ll{R; l"^ (Z)) . 
Observe that the vectors C^^Pj, j = I, . . . , N are the columns of the pre-Gramian 

(2.14) J$ = . . . 

The map Ch links fo-shift-invariant subspaces of L^(M) with /i-doubly-invariant 
subspaces of L^(M;£^(Z)). We recall that a subspace of Lf^{R; i"^ (Z)) is /i-doubly- 
invariant if it is invariant under pointwise multiplication by e^'^**^^ , fc e Z. Obviously 
a subspace S of L^(]R) is to-shift-invariant if and only if the space 

(2.15) Ch{S) = {C,JJ eS} 

is /i-doubly- invariant. T.P. Srinivanasan gave a characterization of doubly-invariant 
spaces in terms of range functions (see jH],[S]). As remarked by de Boor, DeVore 
and Ron [BDR] a similar characterization of shift-invariant spaces follows from it 
(see Proposition |2.1| below) . In our context a range function is a /i-periodic map 
TZ from K to the closed subspaces of £'^{Z). The function TZ is measurable if the 
map V which maps a point a; G K to the orthogonal projection 'P{x) onto TZ{x) is 
weakly measurable in the operator sense, i.e. the function 

X (P(a;)(^,'(/')f2(z) 
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is measurable for all ip and -i/" G Range functions which are equal almost 

everywhere are identified. 

Proposition 2.1. A closed subspace S of L'^(M.) is to-shift-invariant if and only if 
(2.16) S={f e L^(R) : C,Jix) e Uix) for a.e. x eR} 

for some measurable h-periodic range function TZ. 

Obviously to-shift-invariant subspaces with the same range function coincide. 



This observation shall be used in the proof of Theorems 3.6 and 3.7 



In Theorem |2.2| below we compute the range function of the space of band- limited 
functions 

B^ = {fe L2(M) : suppif) C [-0,0]}. 

In general it is not a simple matter to compute the range function of a space. 
Fortunately, if the space S is finitely generated, the range function can be written 
in terms of the generators; indeed in |BDR| it was shown that if 5 = S^^to Etnd 
$ = j = l,...iV} then 

(2.17) TZh{S^^tJ{x)^sp&n{C,,^{x):j^l,...,N} 

is the space generated by the columns of the pre-Gramian matrix ^ . This result 
holds also for countable sets of generators. Of course a shift-invariant space S can 
have more than one family of generators; the smallest number of generators is called 
the length of the space 

(2.18) len,^(S') = min{tt$,5 = S^,tA- 

In |BDR] it has been proved that if S is finitely generated then 

(2.19) leUj^ (5) = ess sup {dim7^/,(5<I>,tJ(a;) : xe[0,h]}. 

To state the next theorem we need some notation. We denote by [a] the greatest 
integer less than a. Recall that h — 2TT/to and set £ ^ [^] -f 1. 

If f < /i< then < -0 + ih < - {£ - l)h < h. We denote by /_, /, /+ 
the intervals defined by 

(2.20) /_ = (0,-0 -f£/i), I = {-0 + eh,0~{e-l)h), 1+ = {0-{£-l)h,h) 
Similarly if < h < then < - (^ - l)h < -0 + £h < h; in this case we 

« 2 

denote by , K, the intervals defined by 
(2.21) 

X_ = (0,0- (£- l)/i), K={0-{l-l)h,-0 + £h), K+ = {-0 + £h,h). 
Recall that {Se : ^ e Z} is the canonical basis of ^^(Z). 

Theorem 2.2. Let to be a positive parameter, h — 2Tr /to and set £ = [^] -|-1. Then 
(i) if J < h < the range function of the space B^^ is 



span{(5j 

(2.22) nh{B^){x) = { span{5j 

span{(5j 



"(^- 1) < j < ^- 1}, xcl^ 
-£<j <£-!}, xel 

-£<j<£~2}, xel+. 



Note that if h ~ j then the intervals /_ and are empty. 
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(ii) If : 
(2.23) 



< h < the range function of the space is 



'Rh{B^){x) 



span{(5j 
span{(5j 
spanj^j 



-(^- 1) < i < ^- 1}, x^K^ 
-(^- 1) < J < ^- 2}, x(^K 
-e<j <e-2}, xeK+. 



Note that if h = then the interval K is empty. 

Proof. For the sake of simplicity we prove the theorem only for ^ = 2, i.e. 0/2 < 
h < 0. The proof in the other cases is similar. Denote by Ai the space of all 
functions g in Ll(R;e'^(Z)) such that g{x) G Tl{B^){x) for a.e. x € [0,h]. Let Q 
denote the orthogonal projection of Lf^{^; ^^(Z)) onto A4. By jHl Lemma p. 58] for 
g e Ll(R;e^(Z)) we have 

(2.24) Qg{x) = Q{x)g{x) a.e. x 

where Q{x) is the orthogonal projection of £^(Z) onto TZ{B^){x). Thus, to determine 
TZ{B^){x) we only need to describe the projection Q{x) for a.e. x. Denote by A0 
the space of functions in L^(M) with support in [—0,0] and by P : L^(M) 1-^ A0 the 
orthogonal projection onto A0, that is P/ — X[-0,0]f- By Proposition 2.1 Chi^^) — 
Ch{Bu>) ~ M. Thus we have 

(2.25) Q = ChVL-^ 

because Lh is an isometry. Let = X^nezV'n^ri G i?i(I^; ^^(^))7 then ^^^^{x) = 
ipnix) for X G [nh, {n + l)h], n G Z; i.e. 



E 



X[nh,{n+l)h]'^n- 



Hence 
(2.26) 



Since | < /i< 



P'^/j V = ^Xl-0,0]^l7ih.Xn+l)h]i^n- 



X[—0,0]X[nh,(n+l)h] 



= < 



X[o,/i] 

X[-h,Q] 

X[h,0] 

X[-0-h] 





if n = 
if n = -l 
if n = l 
if n = -2 
otherwise. 



It follows that 

(2.27) P^h V = X[-0 -h] ^-2 + Xl-h,o] V'-i + X[o,h]V'o + X[/i,0] V'l 

Now we find ChPCj^^ip{x)- 
First suppose that h > ^0; then < 



follows that for each -0 G £ (Z) 



■ h < 2h - < h. Hence from (2.27) it 
V'-i(a;) (5_i + ■00(2;) ^0 + "01(2;) (5i a; G (O,0-/i) 



/:,,P£^V(a;) - < ^-i{x) 5-1 + Mx) So 



^ip-2{x) 6-2 + i'-lix) 5-1 + ^o{x) So 



X e {0- h,2h- 
X e (2/1 - 0, h). 



3) 



8 



VINCENZA DEL PRETE 



Suppose now h < |0; in this case 2h — 0<0 — h<h. Hence from (2.271 we have 



ip-i{x) (5-1 + ijja{x) 6o + ^i{x) 6i x e (0, 2h - 0) 



X e {2h - 0,0 - h) 



_^'4!-2{x) 5^2 + i^-iix) 5-1 + ipo{x) 5o xe {0 — h,h). 

Note that if h = 0/2 the intervals {0,2h — 0) and (0 — h,h) are empty. This 
determines completely the projection Q{x) by (2.241 and ( 2.25[ ). □ 

The following corollary is a straightforward consequence of Theorem |2.2[ 

Corollary 2.3. Let I — [^] +1. Then the length of B^^, as t^, -shift-invariant space, 
is 



(2.28) 



2£ if $<h<j^, 



2£-l if 



< h < 



Proof. The thesis follows immediately from (2.191 and Theorem 2.2 



□ 



3. Frames for the space 

Let $ be a finite family of generators for B^^ . In this section we find conditions 
under which is a frame for B^^. First we prove a representation formula for 

the frame operator analogous to the formula proved by Heil and Walnut for Gabor 
frames [HWi Theorem 4.2.1]. From this formula we deduce a simple necessary 
condition (see Proposition 3.2 1. The formula will also be useful in Section |4] to find 
the dual generators. 

Let <I> = {</5j, 1 ^ J < N} be a family of functions in B^, to a positive parameter 
and h ~ 2tt /to. Denote by Q.^ ^ the function 



N 



(3.1) 



r!5 



fc e z. 



Theorem 3.1. // Efc ll^ltJU < 00 then the operator T^ ^^ from f(Z;C^) to 
Bit, is bounded and 



(3.2) 



Tkh 



on L'^(IR), where the series converges in the operator norm. 

Proof. First we show that if / G L'^{M) then [/, LpJ\ is in Ll for l<j<N. Indeed 
\[l^]\^dx^h f + ^h) ^(x + £h)) [I ^] {x)dx 



(l+l)h_ 



f{z) Lpj{z) [f,LPj]{z) dz 



[Y1 ^(z) lp{z + kh) f{z + kh)f{z) dz. 
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Note that we may exchange the sum and the integral because / ipj has compact 
support and the sum is finite. By summing over j and exchanging the sums we 
obtain 



Therefore, by (3.1 1 
(3.3) 



Hence 



By Schwarz's inequality the right hand side is less then II/II2 X^A; 11^* t 
[f,<P]] G L\, for 1 < i < N. By (2.7 1 the Fourier coefficients of [/,^] are 
{f,Ti^hfj)i k E Z. Hence, by (3.3 1, Plancherel's formula and (2.3 1 

AT 

(Ef^l.t„Tfc,j,7) = EEi(/'^fc''^^)i' 

k j — 1 k 



This proves (3.2 1. Hence T^ ^^ is bounded from to e^{Z;C'^) and 
bounded from~F(Z;C^) to bI. 



□ 



Now by using formula (3.2) we show that the operator ^ t unitarily 

equivalent to the operator of multiplication by the matrix G^,t^ acting on the space 
f^(Z)). Our proof, given for the sake of completeness, is a simple alternative 
derivation of a result of Ron and Shen for general shift-invariant spaces \QS^ . 
Let / be a function in L'^(R). By for each j e Z 



J'T^^.T^^.^J'-'fix + jh) = + J'')fix + + kh) 

k 

= Y,ni-l{x + jh)f{x + ih). 



Now we observe that 



ni:t{x+jh)^%-i{x+eh) 



Therefore by (2.11 1 and (3.1 1 we get 

TT^^,Tl,^T-'f{x + jh) - ^(G 



for a.c. X eR j,£ e 7L. 
Jx)fix + eh) j,£eZ. 



By (|2.13| this formula implies that 



(3.4) 



[0,h]. 



This shows that the operator ^ ^ is unitarily equivalent to the operator G^^ta 
defined by 



(3.5) 



S<S',t,gix) = G$^t„(a;)g(a;) 
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for almost every x G [0, ft.]. An operator of this form is said to be decomposable 
into the measurable field x i— > G<i>,ta{^) of operators on We shall use this 

representation in Section 3 to find the dual generators of frames. 

In the following proposition we give a necessary condition for to be a frame. 

The proof mimics closely an argument of Heil and Walnut for Gabor frames [HWj . 



Proposition 3.2. Let $ = {(pj,l < j < N} be a family of functions of 
Eis>,ta ^•s frame for B^^ then there exist < (5 < 7 < oo such that 



If 



N 



(3.6) 



< 



a.e. m 



Proof. Observe that ft-X^^Li — ^% ta- We only prove that < essinf fi^ ^^ be- 
cause the proof of the other inequahty is analogous. Suppose that ess inf fi^ j =0; 
then for each e there exists a set C (0,0), of positive measure, such that 
^% t (^) < ^ f*-"^ ^-^^ X G E^. We may suppose that there exists an an interval 
/ of measure h such that E^ C /. Let / e B^ be defined by / = ; by the 
Parseval and the Plancherel formula 



I (/, Tkt 



E 

Therefore 

N 



.^.)i^ = E 

k 



f ^ 



^2TTik^ 



'» dx 



^ h 



dx 



j = l,...,iV. 



EEk/'-^*»^.)|'- / 

.7 = 1 k "'^ 



XEAx)nl,{x) dx < eh\\x,e ^ ehWfW 



This contraddicts the fact that E^ t is a frame 



□ 



Next, we give a necessary and sufficient conditions for E^^ta to be a frame or a 
Riesz basis for B^, when $ is a subset of B^^ of cardinality len^ (B^) (see Theorems 



3.6 and 3.7 below). Our characterization will be given in terms of the pre-Gramian. 
Strictly speaking, the pre-Gramian J^^to is an infinite matrix. However we shall 
see that all but a finite number of the rows of J^^to vanish. Hence we may identify 
it with a finite matrix. We shall need the following 

Lemma 3.3. Let h be a positive number, £ = [^] +1 and g G B^^. Let /, and 

K^, K, be the intervals defined in the previous section (see {2.20) and {2.21)). 
(i) Suppose that j <h < . Then Tjh g{x) — 



(3.7) 



if xe I- and j ^ {-{£ - I) < j < £ 
if X el and j ^ {-£ < j < t - 1} 
if X e L+ and j ^ {-£ <j<i-2}. 



-1} 



Note that if h ^ j then /_ and /+ are empty 
(ii) Suppose that < h < 



J. Then Tjh g{x) = 
and j ^ {-{£- I) <j <£ 



(3.8) 

Note that if h = ^ 



if X e K- 

if x&K and j i {-(£ -\)<i<l 
if xe K+ and j ^ {-£ <j<£-2}. 
^ then K is empty. 



1} 
2} 



FRAMES FOR BAND LIMITED FUNCTIONS 



11 



We omit the proof which is straightforward. 



We consider separately the two cases h < 



and 



< h. Assume first that 



J < h < j^. Then \en^^{B^) = 2£ by Corollary Let <^ ^ {ipj : I < j < 21) 
be a subset of B^^ of cardinality 2t. By Lemma ta. 3| all the rows of the matrix 
J vanish except possibly {Tjhipi,Tjh'f2, ■ ■ ■ , Tjii(p2e) , —£ l£ j < ^ ^ 1- Thus we 
identify the infinite matrices ^ t t ^'^^ ^* to with their 21 x 21 submatrices 
corresponding to these rows. The entries of the Gramian matrix are 

J^., ^Wk,^] 1<J<2£ l<k< 2L 

By Lemma 3.3 the i-th column of ^ , 1 < i < 2£ is 
(3.9)^ 







T(t-2)hft 

We note that 



in /_ 



•-Pi 



in / 



Vh 



T(i-2)h^i 





in /_L 



(3.10) rankG^ t^ = rank j^. 

We shall use the following result of Bownik |Bl Thm 2.3] which characterizes the 
system of translates E^^t^ ^ being a frame or a Riesz family (for the space it 
generates) in terms of the "fibers" {Ch'fix) : ip € <&}. 

Theorem 3.4. Suppose <& C L^(M") is countable and let H be the subspace of 
i^(M) generated by -E$,t„. Then 

(i) -E$.f„ is a frame for H with constants A, B if and only if {Ch<p{x) : ip G <!>} 
is a frame for 'R-h{H){x) with constants A, B for a.e. x € [0, h]. 

(ii) is a Riesz basis for H with constants A,B if and only if{ChP> '■ p> G *&} 
is a Riesz basis for TZh{H){x) with constants A, B for a.e. x € [0, h]. 

To apply Bownik's theorem in our context we need a simple lemma of linear 
algebra. Let J be a n x m matrix with complex entries, n < m; we shall denote by 
II J|| the norm of J as linear operator from C™ to C" and by [J]„ the sum of the 
squares of the absolute values of the minors of order n of J. 

Lemma 3.5. Let vi, . . . ,Vjn be m vectors in C", m, > n, and denote by J the 
matrix {vi, . . . ,Vm) whose j-th column is the vector Vj. 

(i) // [J]n > then {vi, . . . is a frame of C" with frame constants A > 
[J]n\\Jf^^^"\ B = \\Jf. Conversely, if {vi, . . . ,Vm} is a frame for 
with constants A and B, then [J]„ > A^ and \\J\\ < B^/'^. 

(ii) If m = n and det J > then {vi, . . . ,v„i} is a Riesz basis o/C" with con- 
stants A — det{J)^ \\J\\^^-^^"^ and B — ||J|p. Conversely, if {vi, . . . ,Vm} 
is a Riesz basis for C" with constants A and B then det( J) > A"/^ and 

\\J\\ < Sl/2. 
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Proof. Let T : C™ C" be the synthesis operator associated to {vi, . . . , Vm}, i-e. 
Tz = J2iLi '^j^j ^ ^ We observe that TT* and JJ* have the same 

eigenvalues Ai < • • • < A„, because the matrix J represents the operator T with 
respect to the canonical bases of C™ and C". Since 

Ai / < TT* < A„ / 

{vi, . . . ,Vm} is a frame for C" if and only if Ai > 0. In such a case Ai is the 
biggest lower frame bound and A„ is the smallest upper frame bound. Moreover 
A„ = ||JJl| = ||Jf. 

Now suppose that [J]„ > 0. By the Cauchy-Binet theorem [J]„ = det(JJ*) = 



^^Aj. Thus 



Ai 



det(JJ*; 

IVk=2 

j} is a frame for C 



> 



\ n— 1 
An 



[J]r 



> 



||J||2(«-1) 

with constants Ai > [J]„|| J|p(i"") and 



,Vm} is a frame for 
Ull < and 



with constants A and 



and {vi, . . 

K = \\jr. 

Conversely, suppose that {wi, . . . 
B. Then A < Ai < A„ < B. Hence 

[J]n = det( JJ*) > A'/ > A" > 0. 

This concludes the proof of part (i) of the Lemma. 

To prove the second part it suffices to observe that {vi, . . . , w„} is a Riesz basis 
of C" if and only if T is an isomorphism and that, in such a case, Riesz constants 
are also frame bounds. Moreover [J]„ = | det Jp when m = n. □ 

Theorem 3.6. Suppose that j <h < j^. Let <i> = {(pj, 1 < j < 21} he a subset 
of B^. Then E^ t^ is a frame for B^ if and only if there exist positive constants 
(5, 7, cr and rj such that 



21 



a. e. m 



(-0,( 



> a 



I det J<j, t I > ry 



a.e. m 



I. 



(3.11) 

(3.12) 

(3.13) 

Ifh = 

Bi^ if and only if conditions (3.11) and (3.13) hold. 

Proof. First we shall prove the theorem for j < h < j^^^- 

Let E^^ta be a frame for B^ with frame constants A and B. This implies in 
particular that B^ coincides with the space S^^ta generated by E^ ta- Condition 
(3.111 follows from Proposition 3.2 Thus we only need to prove (3.12) and (3.131. 
We recall that the columns of the pre-Gramian J^ ^ 



the intervals /_ and are empty. In this case E^_ta a Riesz basis for 



are the vectors Ch^pj, 
j = l,...,2i' by (2.141. Thus, by Theorem 3.4 'i), the columns of J^ ^ (x) are 
a fr ame with constants A,B for the space TlfjB{uj)){x) for a.e. x. By Theorem 
2.2 we may identify canonically TZh{B{uj)){x) with C^^~ "^ for a.e. x £ I- LI 1+ and 
with C^^ for a.e. x in /. Thus, by applying Lemma 3.5 'i) with Vj — Cyjfjix) and 
J = I {x), we obtain that [J^ ^ {x)]2i-i > A^^^^or a.e. a; in /_ U /+ and 

det( J$ fix)) — [>^$ t {x)]^^^ > A^ for a.e. x in /. This proves that conditions 
(3.11 1- (3. 13 1 are necessary. 
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To prove sufficiency assume that conditions (3.11 1-(3.13 1 are satisfied. First we 
prove that the space S^^to spanned by is B^^. Since both are to-shift invariant 

spaces it is enough to show that their range functions coincide almost everywhere. 
We recall that the range Tlh{S'S>.ta) of S'*,t„ is the space spanned by the columns of 

|j I < ^- 1 } and the latter 



J$ J . If a; is in /_ then by (3.9 1 7?.ft,(S'$_t^)(a;) C span{(5. 



space coincides with TZh {B^}{x) by ( 2.22 1 . On the other hand, rank J, 



2£-l 

by (3.9l and assumption (3.12). Thus TZh{S<s,^ta){x) — T^h{Buj){x) because both 



have dimension 2£ — 1 . Similar arguments show that the range functions coincide 
almost everywhere also in /-|_ and in /. 

Next we observe that || ^ (a;)|| < for a.e. x in [0,/i] by (3.111. Moreover 

[J<i, jjx)]2f-i > for a .e. x in /_ U/+ by (|3.12| and [J^^t^{x)]2i ^ | det(J$ ^Jp >0 
for a.e. a; in / by ( |3.13[ ). Thus, by Lemma [3.5[ i), the family {Ch^{x)^ . . . , Ch02t.{x)} 
is a frame for TZh(B^)(x) for a.e. x in [0,h]. The upper frame constant B — 
II J$ ( {x)\\^ is bounded from above by 2ij almost everywhere in [0,h]. The lower 
frame constant A is bounded from below by 

[■h.tAx)hi-i II J$ tA^W^'-^'^ > a(2^7)^'"''^ a.e. in /_ U /+ 



and by 



\detJ^,^ix)\^\\J^,Sx)\\ 



Thus E^^ta is a frame for B^^ by Theorem 
theorem when j < h < 
To prove that \i h = 



2(1-2^) > ^2(2^^) (l-2£) a.e. in/. 

3.4(i). This concludes the proof of the 



£-1/2- 

^ then E^^ta is a Riesz basis for B^ if and only if conditions 

and Lemma 



(3.11 1 and (3.13 1 hold, one argues in a similar way using Theorem 3.4 



3.5 'ii). We omit the details. 



□ 



Remark If £ = 1 condition (|3.12| is superfluous. Indeed, if /i = the intervals 
- and /+ are empty. If ! 



<h <20 then (|3.12| follows from (|3.11|) because 

in /_ 



tJ2£~l 



2 2 



1^1'= 



1^1 



J-hY.]=i\'?3? in/+, 



and the conclusion follows because /_ C ( 
Next we consider the case < h < ^ 

2.3 



-0,0) and /+ C r,i(-0,0). 
J. Then len^ = 2^- 1 by Corollary 

- 2 ° 

Let <& — {ipj : 1 < j < 2^—1} be a subset of i?„ of cardinality 2£— 1. By Lemma 



3.3 all the rows of the matrix ^ , except possibly \Tjh(pi,Tjh(p2T ■ ■ ■ iTjh^2tj 
J ^ ^ ~ 1 vanish. Thus we identify the infinite matrices j t t and G 



with their 2£ — 1 x 2£ — 1 submatrices corresponding to these rows. The i-th column 
of J, 



1 < i < 2^ - 1 is 



(3.14) 



r-{e-l)hfi 



T-{l-l)h'fl 



T(e-2)h'Pi 




iiiK 



T{e.-2)hfi 



inifi 
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Theorem 3.7. Let ^ C such that $ = {ipj, I < j < 21 - 1}, and £ ^ I such 
that <h< -j^. Then E^^t^ is a frame for B^^ if and only if there exist positive 

constants S, 7, a and r/ such that 

2e-i 

(3.15) ^^Yl l^jl^ - ^ (-^^'0)' 

j=i 



(3.16) 
(3.17) 



21-2 



> a 



|det J^ , \ > T] 



a.e. in K, 



in U 



If h = then ta is a Riesz basis if and only if conditions {3.15) and (3.11) 
^2 ' ° L 

hold. 

The proof is similar to that of Theorem |3.6[ We omit the details. 



4. The dual generators 

Let N ~ len^ {B^) be the length of B^ as io-shift-invariant space and let <i> = 
. . . , (/3jv} be a subset of B^^. In this section we shall find the dual generators 
$* when -E$,t„ is a Riesz basis or a frame of B^. With a slight abuse of notation 
in this section we shall denote by <I> the vector ((pi, . . . , (pjv) and by the vector 

It is well known that if E^^t^ is a Riesz basis for S'$,t„ then the Gramian matrix 
is invertible and the Fourier transform of the dual generators are given by 



(4.1) 
where 



= G 



denotes the transpose of the vector v. In Theorems 



4.1 



4.7 



we give 

explicit formulas for the Fourier transforms of the dual generators when E^ to is a 
frame satisfying the hypothesis of Theorems |3.6| or |3.7| The proof is based o n the 
dual Gramian matrix G^^ta representation of the operator ^ ^ . From (3.4) 
we obtain 



1, 



,iV. 



By (2.4) the left hand side is equal to Ch'^k- Hence 



...N 



which, by (2.14), can be written 



(4.2) 



As in Section 2 we identify the infinite matrices and J^. 



with 



N X N matrices by neglecting their vanishing rows and columns (see the discussion 



after Lemma 3.3). Thus we shall interpret (4.2) as an identity between N x N 
matrices. 



Under the assumptions of Theorems 3.6 and 3.7 the interval [0, h] is the disjoint 
union of three intervals where the pre-Gramian is either invertible or has rank A^— 1. 
In the latter case either the first or the last row of the pre-Gramian vanishes. In 
this case we shall denote by J$,t„ the (A^ — 1) x submatrix of ^ obtained by 
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deleting the vanishing row from j . It is straightforward to see that in this case 
equation (4.2 1 reduces to 

(4.3) 



We regard (4.2 1 and (4.3 1 as equations for the unknowns J^, ^ and J$*,to> respec- 
tively. In the intervals where the matrix ^ is invertible we can solve for J^, ^ 



in ( |4.2[ ), obtaining that 
(4.4) 



In the intervals where the rank of ^ is iV — 1 we can solve for the submatrix 



^ obtaining that 
(4.5) 



We recall that if A is a x [N — 1) matrix of rank — 1 then its Moore-Penrose 
inverse IS 



At = {A*Ay^A* 



(4.6) 

(see |BlG] ). Therefore by (|45| 
(4.7) 

We refer the reader to |BIG| for the definition and the properties of the Moore- 
Penrose inverse of a matrix. 

By using (4.4) and (4.5 1 we shall obtain explicit formulas for the Fourier trans- 
forms of the dual generators. For the sake of clarity first we state and prove the 
result for N = 2,3,4. By Corollary 2^ these cases correspond to < /i < 20, 
2 < h < |0 and ^0 < h < respectively. 

Theorem 4.1. Assume that < h < 20 and let <I> denote the vector {ipi, ip2) where 
<^i, </'2 G Slj- If ( 3.11) and (3.13) hold with £ — I, i.e. if E^ t^ is a frame for B^^, 
then 

—± 

DTh^ in [-0,0 - h] 

/i-i||(I)||-2 (j) in {0 - h,h - 0) 

—Dr^h ^ in [h — 0,0] 

{(p2,—^)- Note that if h = the central 



(4.8) 



where D ^ (det J|, jj^^ and 
interval is empty. 

Proof. Assume first that < h < 20. We recall that /_ = (0, ft- — 0), / = (/i — 0, 0) 
and = (0, /i). The pre-Gramian J,j ^ is 



\/ft 







in /_ 



Vh 



Vh 







ill I4 



in /. 



The same formulas hold for J^, ^ with (pi replaced by ip*,j = 1, 2. Therefore 
~ ^ \\/h<^* in /_ 



(4.9) 



/ft (f> in /_ 

/ftr_,,$ in /+ 



VhT^h<^* in /_ 



+ • 
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By assumptions (3.11 1 and (3.13) j has rank 1 in /_U/+ and rank 2 in /. Hence 

J$..t = (J* t V in I- U /+ and 



(4.10) 



Jib* 



J* 



in /. 



First wc find J<[,. ^ in /_ U /+. By using (|4.5| we get 



/ftll*lp 



V^l|r„h*|| 



in /_ 



in Is- 



By (|4.9l we obtain that <i>* = in /_ and r_/,<i>* = / in I,. Since 



T^)J+ = (0 — 0) we finally get 



/l||ci>|| 



h,h- 0). 



Next we find the dual generators in the remaining intervals. By (4.101 
J$. , = (detJ^ i 



V2 - </5]_ 



a.e. in (/i — 0,0). 

By translating and reminding that the pre-Gramian matrix is /i-periodic 

— ± 

Z) T,, ( ^, - ^ ) = D r,, $ 



in (-0,0 -/i) 
in (/i-0,0). 



This completes the proof of the theorem when < h < 20. 

If ft. = one argues as before; the only difference is that now the interval (h — 
0,0 — h) is empty. □ 

To find an explicit expression of the dual generators when iV > 2 we need 
formulas for the rows of the Moore-Penrose inverse of a x {N — 1) matrix of full 
rank. Given a (iV — l)-ple of vectors {Ui, . . . , Um-i) in their cross product is 



N-l 



ei 62 

ul ul 



ttN 



X [/j = C/i X C/2 . . . X Un-i = det 
1 

where the {ej : j — 1, . . . , N} is the canonical basis of . Notice that if = 2 
then X [/ = [/-'-. Given a vector W in C''^ and an integer fee {1, 2, . . . , — 1}, 



JV-l 



we shall denote by X [/^ ^ C/^^ 4— ) the cross product of the {N — l)-ple 
(Ul, Uk-i,W, C/fe+1, . . . , Un~i), i.e. 



N-l 



X U,{Uk^W) = Ui X U2... X Uk-i X W X Uk+i ... X Un-1. 
j=i 



FRAMES FOR BAND LIMITED FUNCTIONS 



17 



Lemma 4.2. Le M he an wx n invertible matrix. Denote by Rj, j — I, . . . ,n, its 

rows and by Cj its columns. Then the columns of are 



(-l)'=+i(detM)~i X Rj 



l<k< 



and the rows are 



1< k<n. 



Proof. Let Mkj be the fcj-cofactor of M. Then the fc-th row of the matrix cof(Af) 
of cofactors of M is 



Lk = (-l)'^+i(Affciei - Mfe2e2 + • • • + (-1)"-U'4„e„) 



= (-iriXi?,. 

i=i 

The first identity follows from the fact that Af^^ = (det M)-^ cof (Af ^). The proof 
of the second identity is similar. □ 

Lemma 4.3. Le A be a n x (n — 1) complex matrix of maximum rank. Denote by 
Aj, j = 1, . . . ,n — 1, the columns of A and by Pk, k ^ 1, . . . ,n ^ 1 the rows of A^ . 
Then 

n-l 

(4.11) Pk = {-l)"[det{A*A)]-^ X Aj{Ak^W) l<k<n~l 

n-l 

where W = X Ai. 

Proof. Since rank A = n — 1 the null space of A* is the space spanned by W. Let 
Af, the matrix obtained by bordering A with the column W, i.e. 

Ab^[ A W ]. 

n-l 

Then Ab is invertible because det Ab = XAj-W=\W\'^. By |BlGl Thm. 8] 



i=i 



a: 



At 



Thus, for every k = 1, . . . , n — 1, the fc-th row of A^ is the fc-th row of Af^^. Hence, 
by Lemma |4 . 2 1 and the anticommutativity of the cross product, we obtain that 



Pfc =(-l)''+i(detAb) ^Ai X ^2 X ... X Ak-i x Ak+i . . . x A^-i x W 



--{-ir\W\'^ X A,{Ak^W). 
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To conclude the proof we observe that \W\'^ is the sum of the squares of the absolute 
values of the minors of order n— 1 of A, i.e. the determinant of A* A, by the Cauchy- 
Binet formula. □ 



Theorem 4.4. Assume that |0 < ft, < and let $ denote the vector {Lpi,(p2,<fi3) 
where ipj S B^^, j = 1,2,3. // assumptions (3.15)-{3.11) hold with t — 2, i.e. 
if i?$^t^is a frame for B^, then the Fourier transform of the dual generators 



(4.12) 



D Th<^ X T2h$ 

E {Th^_y. $) y_Th^ 

D Th% X T_,j$ 

-E T^h% X (T_h$ X $) 

D r_2/i$ X T_,,$ 



in (-0,0 - 2h) 
in (0 ~2h,h- 0) 
in (ft — 0, — ft) 
in (0 - ft, 2ft - 0) 
in (2ft - 0,0) 



where E ~ h (det 



r 



id D = \/ft(det J|, t )-i. iVo^e <fta< if h 



_ 20 
3 



ffte intervals (0 — ft, 2ft — 0) (0 — 2ft, ft — 0) are empty. 



20 

ft, 2ft 



Proof. Assume first that 4^ < ft < 0. We recall that K^,K and iir+ denote the 



intervals (0,0 — ft). 
By (3.141 the matrix , is 



0) and (2ft — 0, ft) respectively, defined in (2.21 1. 



(4.13) Vh 



in 



Vh 







in K 



/ft 



in 



The same formulas hold for J, 
(4.14) J$,t„ = Vft 



^. J with <i> replaced by <&*. Therefore 
and J$._t^ — Vh 



in K. 



By assumptions (3.161 and (3.171 the matrix j has rank 3 in K^LlK-^- and rank 
2 in iC. Hence J^*^ta = ^$,t„^^ in U i4r+ and J^.^t^ = (J5>,i„)^ in 



4.2 



First we find J^, ^ in . By Lemma 
T_,i<I)* = D $ X r/i$, <!)*=£) r;i$ x T_/i<f> 



Dr-h^ X cf> a.e. in ii'^ 



where D = \/ft(det J$ ^ ) ^. By translating the first and the last identities and 
reminding that the pre-Gramian is ft-periodic, we obtain 



(4.15) 



D T/,$ X T2h^ 
D T_2/i$ X T^h^ 



in (-ft,0 - 2ft) 
in (0,0 - ft) 
in (ft,0). 



The same calculation in gives 

T-2h^* = D T_/i<i> X $ T_,i<i>* = I? <i> X r_ 



2h 
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By translating the first two identities we obtain 

(D Th^ X T2h$_ 
D Tft$ X T_/,$ 
D T_2/»$ X T_,j$ 



Next we find J^, j in K. By (4.14) and Lemma 



in (-0, -h) 
in {h -~ 0, 0) 
in (2/i ~ 0,h). 

the rows of Moore-Penrose 



inverse of JJ, are 



4.3 



"VhE (T_/i<i> X $) X $ 
where E = h (dot J$,to J|,.t^)^^- Hence, by (4.141, 



VhE T_/i<i> X (T_;i<i> X <i>) 



-E (r_;,<f> X <f>) X (f) 



([)* 



-i; T_ft$ X (r_,,$ X $). 



By translating the first identity and using the anticommutativity of the cross pro- 
duct, we obtain 



(4.17) 



E X (cf> X Th^) in (0 ~2h,h- 0) 
E (T_/i<i> X $) X T_/i(f) in (0 - h,2h ~ 0). 



The conclusion follows from formulas (4.151, ( 4.16[ ) and (4.171. This completes the 
proof of the theorem when ^ < h < 0. 

li h = ^ one argues as before; the only difference is that now the interval 
(0 — h,2h — 0) is empty. □ 

Theorem 4.5. Assume that | < /i < |0 an d let $ denote the vector {(pi, (p2, fs, V4) 
where ipj G B^, j = 1, . . . , 4. If (3.11)-{3.13) hold with t = 2, i.e. if E^ta is a frame 
for B^, then the Fourier transform of the dual generators $* — {ipl, (p2, (^s-j, fX) is 



(4.18) 



E T_/i$ X (r_,,$ X $ X r,i$) X Th^ 

D r_2/i$ X T^h^ X T/i$ 

E T_2h$ X T_,,$^x (r_2/i$ X r_ft$ x $) 

D r_3,i$ X r_2/i$ X T^h^ 



in {O,2h~0) 
in {2h — 0,0 — h) 
in (0 — h, 3h — 0) 
in (3ft. — 0, 0) 



where E — h (det J$_j^ ^ )^ anrf 13 = h{detJ^ ^ ) ■ The expression o/<i>* in 
(— LjjO) is obtained by reflecting each interval in (4.18) around zero and replacing 
Tjh with T-jh in the expression of in the reflected interval. Note that if h = ^ 
then the first and third intervals are empty. 

Proof. Assume first that | < /i < ^ . We recall that by /_ , / and /+ we denote the 



intervals (0,2ft, - 0), (2ft 
columns of j are 



■ ■ 

^ m i_ 



ft) and (0 - ft, ft) defined in (2.20). By (3.9) the 



in / 



r_2/i5 




in I4- 
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The same formulas hold for J^, ^ with $ replaced by <i>*. Therefore 



(4.19) 



(4.20) 



and J, 



and 



([)* 



in /_ 



in 



By assumptions (3.121 and (3.13) the matrix ^ has rank 3 in /_ U and rank 

4 in /. Hence J$.^i^ = (J|, j^)^ in /_ U /+ and J$.,t^ = "^^.to""^ 
First we find J^, ^ in /; by Lemma 



4.2 



D T_,i (f> X (f> X T/i (f> 

D r_2/i$ X $ X r,,$ 



Th<^* = - D T_2h $ X T_,, $ X $ 



in {2h — 0,0 — h); here D = h (det ^ ) ^. By translating and reminding that the 
matrix j is /i-periodic we obtain 



(4.21) 



L» Tft a)^X T2ft 0? X Tg/j $ in (-0,0-3/1) 

- T_,i $ X Tft $ X T2h $ in (/l - 0, - 2/l) 

D r_2/i $ X T_,i $ X r,, $ in (2/i -0,0-h) 

-D r_3/j $ X T_2h $ X T_,j $ in (3/i-0,0). 



Notice that if [a, fe] is any of the intervals in the r.h.s. of ( |4.21[ ) the expression of 
<I>* in [a,b] can be obtained from that in [—5, —a], by replacing h by —h in the 
translations Tjh, \j\ < 3. 

Next we find the dual generators in the remaining intervals. First let us consider 
the interval /_ . Here by Lemma 
are 

VhE W X $ X Th^, 



4.3 



the rows of the Moore-Penrose inverse of , 



'hE T_h<^ X W X Th<^, 



'hE T_h $ X $ X T¥ 



where W = T_/i<i> X $ X T/i<i> and E = (det I<s>,t„F^ tj ^ ■ using (4.191 and 
translating we obtain 

E ThW X X T2h$ in {~h,h-0) 

(4.22) ^ ^ Ie T^h^ X W X Th% in (0, 2h - 0) 

E T_2h $ X T_,, $ X T_,,iy in (/i,3/i-0). 

Finally we consider the interval /+. Here the rows of the Moore-Penrose inverse of 
JIcj) I are 



'hE Wo X T_h^ X 



'hET_2h<^ XWoX 



'hET_2h^ X T_h^ X Wo 
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where Wq = t_2/i<& X r^h^ X By using (4.20) and translating we obtain 

L (0-3/1, -ft) 



E T2hWo X T,,$ X T2h^ 

E r_/i $ X ThWo X Tft $ 

E T-2h^ X T-h^ X Wo 

Since ThWo = we obtain 

(E ThW X r,, $ X T2h $ 
E T^h'^ X W X Th^ 
E T_2/>$ X T_/i$ X T^hW 



m 
in 



m 
in 



(0-2ft,O) 
(0 - ft, ft). 

(0- 3ft, -ft) 
(0 - 2ft, 0) 
(0-ft, ft). 



By comparing formulas (4.22) and (4.23) we see that the expressions of $* in inter- 



vals symmetric with respect to zero can be obtained from each other by replacing 
Tjh with T^jh (note that replacing Tjh by r^jh changes also the sign of W). 

Formulas ( |4.21[ ), (4.22) and (4.23) give the dual generators. This completes the 



proof of the theorem when | < ft < ^ 



If ft 



one argues as before; the only difference is that now the intervals 



2ft, 2ft — 0), (—ft, ft — 0) and (0 — ft, ft) are empty. 



□ 



Theorems |4.6| and |4.7| below generalize Theorems |4.4| and |4.5| respectively. We 
omit the proofs, which are analogous to the proofs of Theorems |4.4| and |4.5[ We 
recall that K-,K, and K+ denote the intervals (0, 0-(£-l)ft), (0-(£-l)ft, -0+^ft), 
and (-0 + ih, ft) respectively. 



^-T l£ h < and denote by $ the vector (ipi, (p2, . . . , '^21-1), 
2 

= 1,...,2£— 1. If assumptions (3.15)-(3.11) hold, i.e. if 



Theorem 4.6. Let 

where ipj € B^^, j 

E^^ta ^•s o, frame for B^, then the Fourier transform of the dual generators <f>* 
{vhV*2.---V*2t~i) 



«-fe-i 



in TkhK_ for -{£-l)<k<e-l 



-k-2 



-E X T,h^{$^T_khW) in TkhK for -{i-l)<k<i-2 



e-2-k 

1-2 



in TkhK+ for -i<k<l-2 



where W ^ X Tjh^, E = h^^-^ (det J$,t J| ^ )-i and D = h^^i (det J|, ^ )~\ 



j=-£+l 

We recall that /_,/,/+ denote the intervals (O,-0-|-^ft), (-0 + £ft, - (^ - l)ft). 



and 



l)ft,ft) (see (2.201) 



Theorem 4.7. Let j < h < and denote by $ the vector (ipi, ip2, . . . , ip2e) where 



iPj G Bi^, j — 1,...,2^. If assumptions (3.11)-(3.13) hold, i.e. if E.j,t^is a frame 
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for B^, then the Fourier transform of the dual generators = {lp\, (/?2) • ■ • 

E X T,h^{^ ^T^khW) in Tkhl- for -(^-l)<fc<£-l 

j=-c-fc+i 



in Tfc/i / for — i < k < 



e-k-2 

E X $ ( $ ^ r_fe,,VKo ) in nj, 1+ for - £ < fc < £ - 2 

3 = -t-k 

e-1 1-2 
where W = X Tjh^, VFo = X t,7,$. i? ^ h^^^'^ (det J$ f^J|, ^ anrf 

D = /i^-i (det J|,tJ~^ 



5. Sampling formulas for the space 

In this section we shall apply the previous results to oversampling formulas for 
the Hilbert transform sampling and the derivative sampling in i?„ . In the derivative 
sampling formula the coefficients are the values of the function and of its derivatives 
/^"''i ^ ^ j l£ K, at the points of a uniform grid on M. It was ffist obtained by D. 
Jagerman and L. Fogel for K = 1 and by Linden and N. M. Abramson for any K 
|JFj [L] |LA| . Successively J. R. Higgins derived the same expansion formulas by 
using the Riesz basis method [Hi] . 

In |SFj D.M.S. Santos and P.J.S.G. Ferreira have obtained a two-channel derivative 
oversampling formula for i?^^ with < w by projecting both the Riesz basis 
generators of the space B^^ and their duals into the space B^^^ . With this technique 
the projected family is a frame; however notice that projecting the dual of a Riesz 
basis does not yield the dual frame. Thus the coefficients of the expansions of a 
function computed with respect to the projected duals are not minimal in least 
square norm. 

Let to be such that < h < 20 and let <f> = ((pi, (^2) be a vector such that i?$,t<, 
is a frame for S^^. Then by (12. 5|) 



(5.1) 



•^i)Tkt^^i 



yfeB^ 



1=1,2 feez 



By using the Plancherel and the inversion formulas we see that the coefficients 



(5.2) 



{f,rkt^V:) = {Mjf){kto) 



J = 1,2 



are the samples of the functions Mjf — T ^(pj T f at the points kto^ fc S Z. For this 



reason (5.1) is called a sampling formula. These formulas are useful in applications 



when one wants to reconstruct a signal from samples taken from two transformed 
version of the signal. For instance one may want to reconstruct / from samples 
of / and /' (derivative sampling) or from samples of / and its Hilbert transform 
Hf = —iT^^signTf (Hilbert transform sampling). Both are particular cases of 
the family of frames generated by the translates of two functions tpi , 1^2 such that 
'Pi = Xl-0,0], '?2 = wx[-0,0], where m is a function in L°°(M). 
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Proposition 5.1. Let m be a function in L°°(M) and let $ = {(pi, (^2) where 
(5-3) 'Pi=X[-0,0] = m.X[_0,0]. 

Suppose that < h < 20. Then E^ t^ is a frame for if and only if there exists 
a positive number r/ such that 

(5.4) |m — T_/im| > 77 a.e.in {h — 0,0). 

The Fourier transforms of the dual generators are 



,, in [-0,0-/1] 



/i(l + |m|2)' 



in (0 - h,h~ 0) (^* 



^ h{7n — T-h fn) ' 

Ifh — then E^_to ^ Riesz basis for 



/i(l + |m|2) ' 



^ h{;m—T-h fn) ' 



m -0, 



in {0 — h, h — 0) 



3.6 



are satis- 



Proof. Since det J|, j = h(fn —r^h^) the assumptions of Theorem 
fied. The expression of the Fourier transforms of the dual generators can be easily 
obtained from Theorem 14.11 □ 



By choosing m{x) = — «sign(x) in (5.3 1 we obtain the Hilbert transform frames 
for B^. For h — the associated sampling formula is known as the Hilbert transform 
sampling formula (see |Hil Ex. 12. 9]). The coefficients of the expansion are the values 
of the function / and its Hilbert transform TLf at the sample points kto, k € I1. 
Denote by sine the function = sin(a;)/a;. 



Corollary 5.2. Let fi,f2 be defined by 
(5-5) ^ = Xh0.0] ^ 



-*X[-0,0]Sign. 



If < h < 20 then -B$,f„ is a tigth frame for B^^. The dual generators are ip* — 
{2h)~^ipi for i = 1,2. If h — then E^^t^ is a Riesz basis for B^^. Moreover for any 
f G B^ the folloujing Hilbert transform sampling formula holds 

X ^ (/(fcio)T_fct„ cos (y ) sinc(y ) - {nf){kto)T^ku sin ( y) sinc(y )). 

fcez 



Proof. The assumptions of Proposition 5.1 are satisfied. Thus i?4>,t„ is a frame 



and the expression of the duals follows immediately. The frame is tigth because 
TT* = ^I, since = ^<i>. 



Standard calculations show that if ipi and ip2 are the functions given by (5.5) then 

72 



(5.6) (^i(a;) = ^0cos (— ) smc( — j 



(P2{x) 



. ,0X. . ,0X. 

<sm( — )smc(— ). 



Moreover the coefficients of the expansion formula (5.11 are \/2'k f{kto) and 



'2Ti{'Hf){kto), k gZ. This proves also the expansion formula. 



□ 



By choosing m{x) = ix in (5.51 we obtain the derivative frame for B^. Given a 
function g we shall denote by gs the function g{6x). Note that Jsinc^ = y/^X[-5 sj ■ 
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Corollary 5.3. Let <^i,</32 be defined by 
(5-7) ^ = Xh0,0] ^ = *2;x[-0,0]■ 

// < ft, < 20 then E^^ta o, frame for B^; if h — then it is a Riesz basis for 
Bi^. The Fourier transforms of the dual generators are 

-i(l_M) h-0<\x\<0 



\x\ < h-0 



sign(x) h-0 <\x\< 
\x\ <h- 0. 



Moreover for any f G the following derivative sampling formula holds 



(5.8) 



where 



f = V2^J2 (f{kto)r-ktM - nkt,)T_kt^^*X 



=e I'l * (/i - 0)^sinc/i_0(x) + -^(0811100 - (ft - 0) sinc;i_0) * sinc^i^ix) 



ipl{x) 



(ft-0)-sinC;,_^(x) + 



V2 1 



'COSh_0)(a;). 



Proof. By Proposition 5.1 i?$,t„ is a frame for all < ft < 20 and it is a Riesz basis if 



h ~ 0. From (5.2) we see that the coefficients of the expansion (5.11 are \/2tt f{kto) 
and —^/2Trf'{kto). The expression of the dual generators can be obtained from ipi* 
and (p2* by computing the inverse Fourier transform. □ 



A simple calculation shows that if ft = then 

tpl{x) 



1 ■ 2/02; 1 . 2/02; 

sine (y) ip2(x) = --7^ xsmc { — ). 



/2tt 



Figure 1 and Figure 2 below show the Fourier transforms of the dual generators in 



the Riesz basis case ft 



1, and in the case 



1 and ft 



respectively. 




-1 ^,5 



■1 -0,8 ^.S -0,4 -0.2 0.2 0,4 0,6 



-1 fl,5 



Figure 1: ip^, ^(p2 for h — — 1 Figure 2: (p^, j(p2 for h — ^0 



3m _ 3 

2 
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Remark In ^JF Jagerman and Fogel proved the following two-channel derivative 
sampling formula in the case h = i.e. = ^ 



/(x) = ^/(fcto)(sinc2^(x-fcto) + -f'ikto) sinUx 



ktn 



s\nc-[x 



ktn 



(see also |Hil p. 135]). Thus formula (5.8 1 is an extension of the case = ^ to all 
values of io G [f , ^) (i-e. for all h such that < h <20). 

Our last example is a three channel derivative oversampling formula. To obtain 

a frame with three generators we must choose |0 < /i < i.e. to £ '^^ 
Corollary 5.4. Let ipi,ip2,^P3 be defined by 



<^3 



If h < then E^j^ is a frame for B^; if h ~ ^ then it is a Riesz basis for 

B^. The Fourier transform of the dual generators are 



(5.9) ^*{x) 



2h 



^3 {x^ + ihx + 2/i2, -~i{2x + 3/i), -1) 



{Ah{x),Bh{x),Ch{x)) 
{A_h{x),B_h{x),C.h{x)) 



2,hx + 2h'^,-i{2x-'ih),-l) 



where 
Ah{x) 



1 [x + h) + {2x + h){x + hf 
h? I + {2x + HY + x'^{x + h)'^ 



Bh{x) = 



-0 < X < - 2h 
0-2h<x<h-0 
h — 0<x<0 — h 
0-h<x<2h-0 
2h - < X < 

-i 1 — x{x + hY 

I + {2x + hy + x'^{x + hf 



Ch{x) 



1 h + 2x + x{x + hY 
l + (2x + /i)2 + a;2(x+/i)2 



Note that if h = |0 the intervals (0 ~ h,2h ~ 0) and {2h ~ 0, h — 0) are empty. 



Proof. The family E^,t„ is a frame for Bi^ by Theorem 4.4 To compute the dual 
generators we used formula (4.12). □ 

Thus for each value of the parameter h in [|0, 0) the family E^ to is a frame and 
every signal in B^^ can be reconstructed from the values f{kto), f'^^\kto)^ f^-^^kto), k e 
Z, by the following three-channel derivative sampling formula 

(5.10) _L^^(_i)-i/(-i)(H„)^*(x-fcO, 

1=1 kei, 



where the dual generators are given by (5.9). 

In the rest of this section we present two numerical experiments of reconstruction 
of a band limited signal by using formula (5.10 1. We have chosen the signal f{x) — 
;^^(sinc(a;/2))2 (see Figure 3); note that the function / is the Fourier transform 

of the function (1 — |a;|) + , therefore 0=1. 

In the first experiment we take h 
second experiment we take h 

the first case it is possible to find the analitic expression of the functions if* in (|5.9| , 



|0 so that -B$.t„ is a Riesz basis; in the 
j^0 SO that i?$,to is a frame. We observe that in 
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Figure 3: Signal f{x) = -^{smc{x/2)f 



while in the second case they must be computed numerically. Indeed for h = ^ 
one obtains 



1 



f*3 









- '■ I sinc^ 


[r) 




18 / 





xsmc' , 
2tt V3 



1 



2V27r 



x'^ sinc^ 



Figures 4 and 5 below show the functions tpi, v:'2i ^3 J ^3 ^ 




-20D -100 



■0.8 -0.6 -0.2 0.2 0.6 O.B 



1 

J 






r 



m 200 



■0.4 -0.2 0.2 0.4 0.8 




■200 -100 



-0.8 -0.4 -0.2 0.2 0.4 0.8 



Figure 4: h 



Figure 5: j(P2^^3'^ ^ ' 



Since in the case h = ^0 the expression of the duals is known, it is possible to 
write explicitly the sampling formula (5.10): 
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-^(f)i:(-ir[/(^) 



/27r 



3n7r ^ 



+ /( 



3n7r> 



0V9 



X 



3n7r 



inn \ ^ 
/ 



/ 



(2) 



' ' 



2 (a 



3n7r ^ 



This formula was first given by Linden in [L and Linden and Abramson in |LAj . 
see also |HS] . 

In the case h = the dual generators have been obtained by computing 
numerically the inverse Fourier transforms of the functions in (5.9 1: Figure 6 and 
Figure 7 below show the functions ipi,(f2, and ip'^, jip^Tfl', h— j^0. 




200 -0.B 



,4 0.B 



Figure 6: "I>* for h = ^( 



Figure 7: cpl, \ip*2, ip*; h = 



Figure 8 and Figure 9 show the error in the cases ft. = |0 and h = ^0, respec- 
tively. Notice that in the first case the order of magnitude of the error is 10~* 
while in the second case is lO^'^. In second case, since the functions ip* were com- 
puted numerically, to compute their values at the points x — kto we used spline 
interpolation. This accounts for the different order of magnitude of the error. 





Figure 8: The error with h - 



Figure 9: The error with h = 
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